Backscattering Between Helical Edge States via Dynamic Nuclear Polarization 
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We show that that the non-equilibrium spin polarization of one dimensional helical edge states at 
the boundary of a two dimensional topological insulator can dynamically induce a polarization of 
nuclei via the hyperfine interaction. When combined with a spatially inhomogeneous Rashba cou- 
pling, the steady state polarization of the nuclei produces backscattering between the topologically 
protected edge states leading to a reduction in the conductance which persists to zero temperature. 
We study these effects in both short and long edges, uncovering deviations from Ohmic transport at 
finite temperature and a current noise spectrum which may hold the fingerprints for experimental 
verification of the backscattering mechanism. 
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Topological insulators have a bulk band gap, but can- 
not be adiabatically connected to the conventional insula- 
tors without closing the energy gap [TH5 because they are 
characterized by a non-trivial topological number. The 
quantum spin Hall (QSH) insulators belong to the class 
of two dimensional time reversal (TR) invariant Z2 topo- 
logical insulators. They support a pair of gapless helical 
edge states, which give rise to unique, experimentally 
confirmed [SHE] electrical properties. Due to their heli- 
cal nature, QSH edge states are protected against elastic 
back-scattering from TR symmetric perturbations: there 
is destructive interference between two time reversed am- 
plitudes, one where spin is rotated in a, say, clockwise 
manner and a second one with anti-clockwise rotation. In 
the interference term, this amounts to a full rotation of a 
spin-1/2, corresponding to a factor of minus one, which 
cancels the contribution of the direct term. Backscatter- 
ing due to TR breaking perturbations [9 12 and from 
inelastic processes [7J IT3HT8] is however allowed. 

The protection of the edge states in the presence of 
TR symmetry can be exemplified by considering the scat- 
tering due to a spatially random Rashba spin-orbit cou- 
pling [T§1 - |2"2"] . In Fourier space, the matrix element for 
this process includes a factor (k + k'), where k and k' 
are the momenta of initial and final states, respectively. 
For a TR invariant Hamiltonian, the dispersion relation 
satisfies = such that energy conservation en- 

forces k' = —k implying a vanishing matrix element for 
backscattering. In the presence of an in-plane magnetic 
field however, there is a Zeeman splitting 2A = g/ie-B 
between opposite spin directions. For a linearized disper- 
sion e(k) = Hvk this gives rise to a non-vanishing prefac- 
tor (k + k') = 2A/hv and yields a finite mean free path 
I oc 1/A 2 . This mechanism can explain the experimen- 
tally observed dependence of the quantum spin Hall edge 
conductance on an in-plane magnetic field [6] . 

In this letter, we discuss the influence of dynamically 
polarized nuclear spins on the conductance of QSH edge 
channels depicted in Fig. [T] We argue that flip-flop scat- 




FIG. 1. Quantum spin Hall edge states dynamically polar- 
ize nuclear spins via a hyperfine interaction. Although the 
edge states are drawn as sharp lines, their wavefunctions may 
overlap a large number of nuclei. Only one edge is shown, but 
there is an equivalent and independent transport channel at 
the opposite edge, where the spin polarizations are reversed. 



tcring between electronic and nuclear spins creates a dy- 
namic nuclear polarization, which has the same effect as 
an external Zeeman field and gives rise to backscatter- 
ing between helical edge states. The magnitude of the 
nuclear spin polarization is determined by the ratio of 
bias voltage and temperature. Thus, for a long edge of 
length L ^> £, the local temperature profile governed by a 
balance between Joule heating and electronic heat trans- 
port determines the conductance. We find a non-linear 
current-voltage characteristic at finite temperature and 
a peculiar relation between applied voltage and noise 
power, which together with hysteresis in the current- 
voltage characteristic are signatures of backscattering 
due to a dynamic nuclear polarization. Importantly, this 
mechanism stays effective in the low-temperature limit, 
similar to spin-flip scattering in v — 2 quantum Hall 
edges [23] . 

At low energies in the presence of a nuclear spin polar- 
ization, the helical edge electrons have an effective Dirac 
spectrum governed by 

H = f dxii{x)a z ap (-ivW x + \A{I)\ ^(x), (1) 
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where ipa^ i{ x ) 1S a fermionic field corresponding to a 
right moving (spin-up) or left moving (spin-down) elec- 
tron with Fermi velocity v, a z is a Pauli matrix in spin 
space, and repeated indices are summed over. The po- 
larization of nuclear spins is generated via the hyperfinc 
dynamic Hamiltonian |24| 



M x i) ( 2 ) 



where I i is the spin operator for a nucleus located at 
position n = (x i ,y l ,z i ), a ± = (a x ±ia y )/2, and p n is 
the effective one dimensional density of nuclei seen by 
the edge electrons. Spin flips are generated by this ef- 
fective contact term, where the strength of the hyperfinc 
interaction, A, can be measured using optical or electron 
paramagnetic resonance experiments |25j . Backscatter- 
ing between right movers and left movers can be induced 
by a spatially dependent Rashba spin-orbit coupling term 
that can arise as a result of electrostatic potential disor- 
der [T5H22) . The interactions are described by 



Hr= I dx^ a {x)a v al3 {a{x) 1 id x \^ l3 (x) 
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where the strength of the inhomogeneous coupling is 
set by the scale of its local correlations: (a(x)a(x')) — 
Vr5(x — x'). The total Hamiltonian is given by H = 
Hq + Hf\i p + iJ K where the effective Zeeman gap orig- 
inating from the polarized nuclei A(I) will need to be 
determined from the steady state of nuclear polarization. 

To this end, we begin by calculating the net scattering 
rate per unit length, T, describing the transfer of elec- 
trons between the right and left moving spin polarized 
edge states, \k f) and | — k' J,), due to the interaction 
between the electron and the nuclear spins -ffflip- Em- 
ploying Fermi's golden rule we obtain [2l)H2"5] 
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(1/2 + </(«)» /_(e)[l-/ + (e)] , 
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where /±(e) are the distribution functions for the right 
and left moving edge states, respectively. 

Because each scattering event is associated with a 
nuclear spin flip, the time evolution of the nuclear 
spin polarization can be described with a rate equa- 
tion d t (I(t)) = T/p n - T d (I(t))/p n , where r d is a phe- 
nomenological nuclear spin relaxation rate encompassing 
all other mechanisms that lead to nuclear spin flips. In 
quantum Hall systems, the nuclear spin relaxation rate is 
typically much smaller than the spin flip scattering rate 
that originates from the dipole interaction between the 
nuclear spins and the spin polarized edge states [23J . Due 
to their similarity with the QSH edge states considered 
here, we expect that <C T. 



The rate equation for the nuclear spin polarization 
can be solved for arbitrary electron distribution func- 
tions /±(e). Assuming that the nuclear spin polariza- 
tion is initially zero, the time evolution is described by 
(/(*)) = (7)(1 - e-*/ T »), where 
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It is difficult to obtain a reliable estimate for the char- 
acteristic time-scale r„, where the steady-state nuclear 
polarization (I) is achieved. Based on the experiments 
in the quantum Hall systems [23] , we expect that the first 
term in t^ 1 dominates and r n is on the order of tens or 
hundreds of seconds. 

To evaluate the effect of a finite nuclear polarization 
(I) on electrical transport, we first concentrate on the 
case of a short edge, where the distance between con- 
tacts L is much shorter than the mean free path associ- 
ated with the spin flip scattering. The distribution func- 
tions are then determined by the Fermi functions f± (e) = 
/o(e - Pl/r,T), where f (e,T) = [1 + exp(e/fc B T)]- 1 , 
and Pl/r are the chemical potentials of the left /right 
reservoirs, respectively. At zero temperature (I) = 
i n _l_ rdSTTfi-D SQ t j iat m absence of relaxation, 

the nuclear spins become fully polarized. For Yd = 0, the 
temperature dependence of the nuclear spin polarization 
is described by (I) = h tanh [(p L — p R )/2k B T] . 

We now turn to evaluating the backscattering rate due 
to the random Rashba spin orbit coupling in the pres- 
ence of a dynamically generated nuclear spin polarization 
(I). Denoting initial and final electronic states by their 
momenta \k f) and i), we find the matrix element 
(k'\Hft\k) = —i(k + k') a(k' — k)/L. Taking into account 
that due to energy conservation (k + k') = —A(I)/vH, 
we find for the disorder averaged scattering rate Wkk' = 
2irV R A 2 (I} 2 5(e k - e k ,)/{h 3 v 2 L). The total Rashba scat- 
tering rate is then 
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The edge conductivity due to Rashba scattering is a — 
(e 2 /Trh)£ with I — vtr/2 (the extra factor 1/2 accounts 
for the difference between single particle and transport 
relaxation time). It is related to the two-terminal con- 
ductance via 



G = 



(6) 



In the limit of a short edge with L<l, the conductance 
is changed by 



SG 



e \ 42/ na V * 



(7) 
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FIG. 2. An outline of sequential scattering processes for the 
helical edge states. It is assumed that the nuclei are unpo- 
larizcd ((J) = 0) in the leads but are dynamically polarized 
((/} 7^ 0) between them in a region with potential disorder 
giving rise to random Rashba scattering. 



This equation can be directly verified by calculating the 
change of current caused by scattering between the spin 
polarized edge states. Here, we assume that the phase 
space for scattering is not changed by the presence of 
a spin splitting due to the nuclear polarization. This 
can be understood by calculating scattering states from 
the time independent Schroedinger equation describing 
energy eigenstates as pictured schematically in Fig. [5] 
We assume that there is no spin splitting in the reservoirs, 
and that it is switched on and off in regions of the edge 
where there is no Rashba scattering. Since the nuclear 
magnetization only couples to the electron spin, which is 
assumed to be a good quantum number, a finite (J) leads 
to a change in wave vector, such that Rashba scattering 
has a finite matrix element. After the scattering, (/) is 
switched off again before entering the reservoir, so the 
initial and final momentum are equal in magnitude. 

For long edges with L £, the distribution functions 
of the electrons are no longer determined by the reservoir 
distributions but have to be calculated from a solution of 
a kinetic equation. In the following, we assume that the 
Coulomb interaction between counter-propagating edge 
states gives rise to an inelastic mean free path approxi- 
mately equal to the elastic one |29J . Thus, it seems rea- 
sonable to assume that a long edge region can be de- 
scribed locally by equilibrium distributions with position 
dependent chemical potential n(x), temperature T(x), 
and nuclear polarization (I)(x). Neglecting the relax- 
ation of nuclear spins, the nuclear polarization is related 
to the difference in chemical potential between left and 
right movers, Sfi, which in turn is determined by the elec- 
trical current j via Sfi — jh/e. As current is conserved 
along the edge, we have 

™*> - ^ anh (ssy ■ <8) 

Expressing the electron distribution function as 

f(k,x) = f (Z k ,T(x))+5f(k,x) , (9) 

with £fc = efc + sign(k) A(I) (x) — [i(x) subject to the 
linearized steady state Boltzmann equation Vkd x fo — 



-Sf2/TR, we find the solution 



6f(k,x) = ^<5(a) 



sign(k)Ad x (I) - d x fi - 7^d x T 



(10) 

The current is obtained as j — e J ^VkSf, giving rise to 
a position dependent conductivity 



a(x) 



e 2 h 4 v 4 



nh 2V R A 2 {I) 2 {x) 



(11) 



The conductivity depends on temperature via the expec- 
tation value of the nuclear polarization in Eq. ([8]). In 
principle, interaction effects give rise to an additional, 
explicit temperature dependence of the backscattering 
rate, which is not included in the following. In order 
to determine the conductance, knowledge of the temper- 
ature profile T(x) along the edge is needed. In the limit 
of small bias, the edge temperature is given by the bath 
temperature, T{x) = T. Linearizing Eq. wty we obtain 







2h 4 v 4 fek B T 



ttH Vr \ AttH 



dx 



(12) 



a nonlinear current voltage relation, where in the limit 
a/L < e 2 /h one finds j cx T 2 / V /3 - 

In general, the temperature distribution T(x) must 
be calculated self-consistently, taking into account Joule 
heating due to the transport current [301 131] ■ Denot- 
ing the heat current by jq, it is related to the lo- 
cal temperature gradient by the heat conductivity k = 
aT (k-Q / e) 2 {it 2 / 3) as j'q = —k8 x T. Joule heating deter- 
mines the divergence of the heat current according to 
9 x Jq = j 2 /cr, such that 



a{x)d x [a(x)d x T 2 {xj\ 



J 
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(13) 



Together with Eqs. (|Sj) and (111 this constitutes a nonlin- 
ear differential equation for the temperature profile. We 



now discuss the solution of Eq. ( 13 ) in the limit where 



the bath temperature can be neglected as compared to 
the temperature generated by Joule heating, e.g. we use 
the boundary conditions T(0) = T{L) = 0. Introducing 
the length scale £ = (hv) 4 j \2VrA 2 ) , we define a dimen- 
sionless coordinate x = x/£q and temperature variable 
t(x) ee T(x)/e(j) with 9(j) = jirh/ek B . Eq. Q may 



now be numerically integrated yielding the dimension- 
less temperature profile and nuclear polarization shown 
in Fig. |3] for fixed L/l > 1. 

The chemical potential difference A/j, between the left 
and right reservoirs is related to the current via 



J 



.ttH 
4e 



dx tanh 



t(x) 



(14) 



Since the functional form of the rescaled temperature 
profile t(x) only depends on the dimensionless param- 



eter L/Iq through the boundary conditions of Eq. (13 1, 
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FIG. 3. The induced dynamic nuclear polarization (I)(x) 
computed from Eq. Q for a long (L/£o S> 1) edge using the 
dimensionless temperature profile shown in the inset. 



FIG. 4. The Fano factor calculated from Eq. |l5| is found to 
be parametrized by the finite size scaling form F = 1.219(2) — 

2.03(3)(VL)°- 65(3) . 



it is clear from Eq. ( 14) that the resistance of a long edge 



is Ohmic in the limit of zero bath temperature. However, 
we obtain an interesting scaling irhG/e 2 oc (Iq/L) - 35 of 
the conductance with the length of the wire. The current 
noise can be obtained from [301 [31] 



S 



dx u(x) k^T{x) 



(15) 



and can be expressed in the form S = F GV, where F 
denotes the Fano factor and G is the conductance ob- 
tained from the current voltage relation Eq. (14 1. By 



studying the behavior of Eq. (151 in the long edge limit, 



as shown in Fig. |4j we extract a numerical result for the 
Fano factor of F = 1.219(2) where the number in the 
bracket indicates the uncertainty in the final digit. By 
comparison, the Fano factor for a one dimensional diffu- 
sive wire much longer than the inelastic mean free path is 
V3/2 [32l E3] , if cooling due to electron-phonon coupling 
can be neglected. 

Before attempting to apply the results discussed here 
to the HgTe quantum wells studied in Ref. [5] it is useful 
to put them in the context of other potential backscat- 
tering mechanisms that have been proposed to explain 
the observed deviation from ballistic edge transport. In 
the presence of both time reversal symmetry breaking 
and disorder, the helical edge states will be localized, as 
studied in the presence of magnetic fields and impurities 
[9], random magnetic fluxes 10J and magnetic impurities 
[Til [HI- Alternatively, back-scattering becomes possi- 
ble in the presence of dephasing processes, which have 
been simulated phcnomcnologically using Biittikcr's vir- 
tual probes [7J [T3] . Budich et al. [H] found that inelastic 
processes due to electron-phonon interactions in the pres- 
ence of spin orbit coupling cannot cause backscattering 
between helical edge states to leading order. Other stud- 
ies have shown that avenues towards backscattering from 
magnetic impurities alone [15] . electron-electron interac- 



tions [Jj)J [T7] or the loss of axial spin symmetry [TH] are 
all suppressed in the low temperature limit. 

Backscattering from dynamic nuclear polarization per- 
sists to zero temperature. In order to obtain an estimate 
for the mean free path £q, we exploit the interpretation of 
A(I) as a local effective Zeeman splitting. In this way, we 
can determine the parameters in the scattering rate by 
comparison with the effects of backscattering due to an 
in-plane magnetic field in the presence of a spatially ran- 
dom Rashba coupling. Examining Ref. [5] we estimate 
that {hv) A /2V R « 0.14 (/LteV) 2 /im, which in combination 
with a hyperfine constant of A s» 60 fieV for interactions 
of electrons with the fully polarized nuclei of mercury |34j 
produces £q « 40 fim. 
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